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Symbols 

C L lift coefficient 

c atmospheric trace species 

c fluctuation of trace species from mean value 

D fuselage diameter, m 

e partial vapor pressure, hPa 

/ function 

GOES Geostationary Operational Environment Satellite 

INS inertial navigation system 

L lift of airplane, which is equivalent to weight for level flight 

M Mach number 

M } local Mach number 

p free-stream static pressure, hPa 

Ap static pressure position error, hPa 

p' local static pressure, hPa 

p s ref static pressure interpolated to airplane sensor height by using tower measurements, hPa 

p t total pressure, hPa 

q c dynamic pressure, hPa 

q c local dynamic pressure, hPa 

R gas constant 

r recovery factor for sensor 

S wing area, m 2 

TAMMS turbulent air motion measurement system 

T m measured air temperature, °C 

T r recovery temperature, °C 

T s corrected static air temperature, °C 

T* static air temperature corrected for Mach number only, °C 

T f static air temperature interpolated to airplane sensor height by using tower 

measurements, °C 

T t total air temperature, °C 

T x air temperature error due to heat-transfer effects, °C 

air temperature measured at 76.2-m tower level, °C 
U a true airspeed, m-s” 1 

u east-west horizontal wind velocity component, m-s -1 

u p east- west airplane velocity component with respect to Earth, m-s” 1 

v north-south horizontal wind velocity component, m-s” 1 

v a air velocity with respect to airplane, m-s -1 

v p north-south airplane velocity component with respect to Earth, m-s” 1 

w vertical wind velocity component, m-s -1 


v 



fluctuation of vertical wind component from mean value, m-s - 
vertical airplane velocity component with respect to Earth, m-s 
distance of static port in front of airplane nose, m 
height of airplane static port, m 

distance of airplane static port below 91.5-m tower level, m 

angle of attack, deg 

angle of sideslip, deg 

ratio of specific heats 

ratio of specific heat for dry air, 1 .4028 

air temperature probe recovery correction 

pitch angle, deg 

airplane heading, deg 



Abstract 


A turbulent air motion measurement system (TAMMS) was integrated onboard 
the Lockheed 188 Electra airplane (designated NASA 429) based at the Wallops 
Flight Facility in support of the NASA role in global tropospheric research. The sys- 
tem provides air motion and turbulence measurements from an airborne platform 
which is capable of sampling tropospheric and planetary boundary-layer conditions . 
TAMMS consists of a gust probe with free -rotating vanes mounted on a 3.7-m epoxy- 
graphite composite nose boom , a high-resolution inertial navigation system (INS), 
and data acquisition system. A variation of the tower flyby method augmented with 
radar tracking was implemented for the calibration of static pressure position error 
and air temperature probe. Additional flight calibration maneuvers were performed 
remote from the tower in homogeneous atmospheric conditions . System hardware and 
instrumentation are described and the calibration procedures discussed . Calibration 
and flight results are presented to illustrate the overall ability of the system to deter- 
mine the three-component ambient wind fields during straight and level flight 
conditions . 


1. Introduction 

Airborne air motion measurement systems have pro- 
gressed from the research and development stage 
(Lenschow 1972; Meissner 1976; Lenschow et al. 1978) 
to a resource that is appropriate for a variety of applica- 
tions in the atmospheric sciences. The accurate determi- 
nation of the ambient three-component wind field from a 
moving platform has proven to be a very challenging 
problem because a relatively small result must be deter- 
mined from the difference between two large values. The 
relative motion of an aircraft with respect to the Earth 
must be subtracted from the velocity of the air with 
respect to the aircraft in order to obtain the ambient wind 
field. As a result, both sets of these relative velocity mea- 
surements require a high level of accuracy to meet the 
measurement requirements for modem applications. 
Advances in the accuracy of inertial sensors and data 
analysis are responsible for the development of these sys- 
tems and for their expanding utility (Axford 1968; 
Telford and Wagner 1974; Telford, Wagner, and Vaziri 
1977; LeMone and Pennell, 1980; Frost, Paige, and 
Nelius 1991). 

Presently, air motion systems operate on many types 
of research aircraft that encompass a broad spectrum of 
flight conditions. These systems basically follow the 
same approach by measuring relative ground motion via 
an inertial navigation system (INS) but differ slightly in 
their measurement approach for relative air motion. Most 
systems have used either a gust probe (a pitot-static tube 
with flow directional vanes) affixed to a nose boom 
(Lenschow 1972; Gilmer, McGavin, and Reinking 1978; 
Greenhut and Gilmer 1985; Ritter, Smith, and Cahoon 
1987; Anon. 1991b), the radome method (Larson, 
Flecher, and Siemers 1980; Brown, Friehe, and 
Lenschow 1983; Scott et al. 1990; Tjemstrom and Friehe 
1991 ), or a probe-mounted differential pressure flow sen- 


sor (Armistead and Webb 1973; Anon. 1991a; Lenschow 
et al. 1978) for making airflow measurements relative to 
the aircraft. 

In support of the NASA role in global tropospheric 
research, the Lockheed 1 88 Electra airplane (designated 
NASA 429) based at the Wallops Flight Facility was out- 
fitted with a turbulent air motion measurement system 
(TAMMS). TAMMS utilizes a gust probe with free- 
rotating vanes affixed to a 3.7-m nose boom and a high- 
resolution INS used in concert with fast-response chemi- 
cal instrumentation to make airborne flux measurements 
of trace species. The latest technology was employed to 
construct a thick-walled epoxy-graphite composite boom 
system with a natural frequency that is sufficient to pre- 
serve measurements with a 10-Hz bandwidth. The 
TAMMS has obtained flux measurements of trace spe- 
cies over the Amazon River Basin as well as tundra and 
wetland regions in North America (Ritter et al. 1990, 
1992, and 1994). 

The accurate measurement of the three-component 
ambient wind from an airborne platform is complex and 
requires an overall calibration of the air motion system. 
The calibration must cover the entire range of flight con- 
ditions encountered during data flight legs. Systematic 
errors in pressure and temperature measurements are 
induced by the enhanced kinetic energy and distortion of 
the airflow local to the aircraft body, which can often 
contaminate measurements acquired in a flight program. 
Calibration flights and in-flight maneuvers provide data 
needed to characterize and compensate for these errors. 
Errors in the velocity of the aircraft with respect to Earth 
are dictated by the quality of the inertial navigation sys- 
tem. The discussion and characterization of inertial sys- 
tem errors, however, are not within the scope of this 
paper. Instead we focus on the characterization and cali- 
bration of the gust probe and the overall TAMMS 



platform employed for the determination of species flux. 
The following topics are discussed: 

1 . Fundamentals of air motion measurement 

2. Instrumentation 

3. An analysis of the fundamental issues of concern 
for airborne air motion measurements 

4. Flight calibration techniques and validation 
results 

2. Background of Air Motion Measurement 

The three-component ambient wind field vector V, 
as related to a local Earth coordinate system, is defined to 
be 

V = iu + jv + k w (1) 

where u and v are the horizontal wind velocity compo- 
nents and w is the vertical wind velocity component. The 
unit vectors i, j, and k represent the three-component 
direction in a right-hand orthogonal coordinate system. 
The ambient wind field V is obtained from an airplane by 
the vector sum of the air velocity with respect to the air- 
plane (V a ) and the airplane velocity with respect to the 
Earth (V p ) as expressed by 

V = V a + V p (2) 

The transformations and actual expressions used to deter- 
mine the ambient horizontal ( u and v) and vertical (w) 
wind components from airplanes into the standard mete- 
orological frame of reference have been derived and dis- 
cussed in detail by Axford (1968) and Lenschow (1972). 
These complex equations, however, can be simplified 
to the following approximate calculations (Lenschow 
1972): 



u = -U a sin(y + (3) + Up 

(3a) 


v = -U a cos(\|/ + |3)+ v p 

(3b) 


w = -U a sin(0 - a) + w p 

(3c) 

where 

V a 

true airspeed 


V 

airplane heading 


p 

sideslip angle 


0 

pitch angle 


a 

angle of attack 


V P' Wp 

three-dimensional airplane velocities related 

to Earth coordinate system (obtained from 


INS) 


A diagram of the coordinate system, airplane attitude 
angles, and reference axes used to calculate air velocity 
components is shown in figure 1. Equations (3) identify 
parameters that make the most significant contributions 
in the calculation of ambient wind components. These 
equations can be used to estimate the required accuracy 
limits of these parameters for making useful air motion 
measurements. Lenschow (1986) discusses and provides 
an example of these limits for airplanes flying at 
100 m-s _1 , which approximates the typical airspeed of 
the TAMMS platform. 

3. Instrumentation 

The Lockheed 188 Electra airplane (based at 
Wallops Flight Facility) shown in figure 2 has been out- 
fitted with TAMMS which consists of three subsystems: 
(1) a gust probe equipped for fast-response meteorologi- 
cal and flow angle measurements, (2) a high-resolution 
INS, and (3) a data acquisition system. (See the appendix 
for airplane specifications.) The fluctuations in the three- 
component ambient wind field measured by TAMMS are 
combined with measurements of a species made from 
fast-response (>10 Hz) chemical instrumentation to make 
airborne eddy correlation flux measurements. The verti- 
cal flux of a species c is given by 

F = v7? (4) 

where 

w = w — w (5) 

c = c-c (6) 

and the bar over the symbol and prime denote the mean 
value and the fluctuation from the mean value, respec- 
tively. Sample periods are typically 10 to 20 min in dura- 
tion to obtain a statistically significant sampling. 

The flight calibration of TAMMS employed the 
tower flyby method to accurately characterize the impact 
of airplane motion on the measurement of free-stream 
static air pressure and static air temperature. This method 
requires an instrumented tower for correlative static 
ground measurements. A C-band radar facility comple- 
mented the calibration procedure by providing continu- 
ous and accurate altitude data during tower flyby 
maneuvers and additional calibration flight maneuvers 
performed remotely from the tower. The airplane system 
and the ground instrumentation associated with the 
TAMMS flight calibration are discussed in sections 3.1, 
3.2, and 3.3. 
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3.1. Nose Boom and Meteorological 

Instrumentation 

The nose boom is a conical composite structure fab- 
ricated from high-modulus graphite preimpregnated with 
an epoxy resin. It is 3.7 m in length, 25.4 cm in diameter 
at the base, and 7.6 cm in diameter at the tip. The wall 
thickness of the cone decreases from 0.61 cm at the base 
to 0.46 cm at the tip. An NACA research- type combina- 
tion pitot-static tube is attached to the tip of the graphite 
boom. This position places the meteorological instru- 
ments and airflow sensors as listed in table I at a distance 
forward of the fuselage which significantly reduces the 
influence from the airflow around the airplane. Con- 
structing a boom system of sufficient length that will 
totally remove the influence of airflow about an airplane 
is not physically possible. Differential and absolute 
pressure transducers for impact and static pressure mea- 
surements, respectively, are housed in the pitot-static 
tube-boom mounting adapter. Fast-response temperature 
and humidity sensors are mounted to the adapter. Flow- 
direction vanes are located on an NACA combination 
pitot-static and flow-direction probe. The overall instru- 
ment configuration is shown in figure 3 and weighs 
approximately 31 kg. The boom is mounted at the base to 
a support frame housed within the radome which accom- 
modates the INS. This support frame is rigidly attached 
through truss structures mounted to the forward fuselage 
bulkhead. A Honeywell Primus 800 weather radar 
mounts above the support frame between the truss struc- 
tures without field-of-view restrictions in the forward 
or zenith directions. The lowest natural frequency of 
the overall INS, support structure, and boom system is 
15.3 Hz, which is sufficiently beyond the desired 
0-10 Hz bandwidth. 

The static and differential pressures, which are criti- 
cal parameters in determining true airspeed U a , are mea- 
sured by Paroscientific Digiquartz pressure transducers, 
models 215-AW-002 and 5006-D-002, respectively. 
These transducers employ a crystalline-quartz resonator 
which measures change in pressure via their pressure- 
dependent oscillation frequency. Busse (1987) provides a 
detailed explanation of transducer mechanics and opera- 
tion. Application accuracy of the transducers is cited 
by the manufacturer to be 0.01 percent of full scale with 
1 x 10 -8 resolution. Both the differential and absolute 
pressure transducers are housed in a thermally insulated 
and thermostatically controlled “cocoon.” Thermal sta- 
bility of the transducers is obtained with Du Pont Kapton 
laminated etched-foil heating elements controlled with 
an inhouse feedback circuit. The thermal environment of 
the pressure transducers is set at a control temperature of 
44°C±0.5°C to maintain temperature stability during 
normal operating conditions. 


Temperature measurements needed to determine U a 
are made with a Rosemount model 102 non-deiced total 
air temperature T ( sensor housing with a fast-response 
platinum sensing element (E102E4AL) (DeLeo and 
Werner 1960; Werner, DeLeo, and Rogal 1961; 
Stickney, Shedlov, and Thompson 1990). The 
E102E4AL element has a nominal 50-Q resistance and a 
20-ms response time. The Rosemount signal conditioner 
was replaced with an inhouse signal conditioner designed 
to have a 0- to 10-V output. Rosemount total temperature 
sensors normally exhibit a measured temperature T m of 
=0.995T r The sensor is mounted on the starboard side of 
the pitot- static tube-boom adapter. A high correlation 
between the yaw angle and total air temperature T { mea- 
surements observed during in-flight calibration maneu- 
vers indicates a shadowing effect from the pitot-static 
tube and/or adapter, which is discussed later. Also, due to 
the support structure for the temperature sensing ele- 
ment, the T t system actually has a two-term exponential 
response (Lenschow 1972; McCarthy 1973). A numeri- 
cal filter was employed to recover the high-frequency 
signal that is lost because of the presence of the support- 
ing structure for the element (Ritter, Smith, and Cahoon 
1987). 

Airflow angles are determined from a pair of free- 
rotating balsa vanes mounted orthogonally to the instru- 
ment head as shown in figure 3. The vanes are oriented to 
measure angle of attack and flank angle, that is, the inci- 
dence angles of airflow relative to the airplane in the ver- 
tical and horizontal planes, respectively. Sideslip is 
assumed to be equal to flank angle for small angles of 
attack. The vane assembly is attached to a shaft which 
rotates in a hollow strut and flange machined from a 
solid piece of stainless steel for mounting. A bronze 
sleeve bearing supports the outboard end of the shaft 
with a synchro repeater attached to the inboard flange 
end. A 16-bit synchro- to-digital converter was used to 
obtain the data. As a result, airflow angles of 0.005° were 
capable of being resolved. With equation (3c), this flow 
angle translates into a sensitivity of 1 cm-s” 1 for the ver- 
tical velocity component. 

Richardson (1971) discusses the design of the vanes 
and provides data on the dynamic properties as employed 
by the TAMMS on the Electra airplane at Wallops. 
Advantages of ffee-rotating vanes include the measure- 
ment of the airstream angles directly; therefore, the need 
is alleviated for accurate knowledge of aerodynamic 
coefficients for data interpretation, and the elimination of 
acceleration sensitivity by mass-balancing each vane 
(Lenschow 1971; Richardson 1971). The natural fre- 
quency of free-rotating balsa vanes was determined by 
Richardson to be approximately one-half times the 
square root of the dynamic pressure measured in pascals, 
which translates to 25 to 50 Hz throughout the range of 
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the Electra flight envelope. The lower detectable limit of 
the length scale for turbulent eddies would then be 
approximately 4 to 8 m with the assumption of an aver- 
age airspeed of 100 m-s -1 . Spectral analyses of flight 
data from these vanes support these findings and indicate 
good response over the 10-Hz bandwidth of interest. 

A Lyman- Alpha hygrometer (Buck 1976) is 
mounted on the boom tip adapter aft of the pitot- static 
tube and flow angle vanes. A slower response EG&G 
137 hygrometer is mounted in a window port on the fuse- 
lage and used to normalize the Lyman-Alpha signal. 

3.2. Inertial Navigation System 

The Litton Model LTN-72RH gyro-stabilized iner- 
tial navigation system is employed on the Electra aircraft 
for TAMMS. This system is the high-precision version 
of the standard LTN-72 model. The RH model has been 
primarily developed for scientific applications with strin- 
gent requirements such as those needed for scientific air 
motion measurements. It provides an increased update 
rate of approximately 24 data frames/s on both a binary 
and binary-coded decimal (BCD) bus in ARINC 561 for- 
mat. The position error drift is approximately 0.5 km/hr 
(personal communication from Charles Robinson of 
Litton Industries, Inc.). The INS is mounted to the base 
of the boom within the Electra radome by way of an 
environmentally controlled housing. 

The airplane attitude angles obtained from the INS 
were electrically fed through a 16-bit synchro-to-digital 
converter to yield an angular resolution of 0.005°. 
Equation (3c) indicates that for an airplane speed of 
100 m-s -1 , the required angular resolution for 0 is 0.06° 
in order to meet a short-term resolution requirement of 
1 cm-s -1 for the vertical velocity fluctuation. The INS 
used by the TAMMS platform, therefore, easily meets 
this requirement. The vertical velocity of the airplane w p 
is derived by integrating the vertical acceleration output 
of the INS and bounding it by the third-order barometric- 
inertial loop algorithm as suggested by Lenschow 
(1986). 

The long-term accuracy of the horizontal velocities 
u p and v p are dictated by INS drift rate. A thorough dis- 
cussion of inertial systems and the errors present in 
the resultant velocity measurements are presented by 
Broxmeyer (1964) and Kay ton and Fried (1969). 
Lenschow (1972) gives a general discussion on the types 
and orders of magnitude of errors associated with inertial 
systems. 

3.3. Aircraft Data Acquisition System 

The TAMMS data acquisition system consists of 
three subsystems: (1) signal conditioning and interface 


electronics, (2) a pulse code modulation (PCM) system, 
and (3) a data recording and display system. A block dia- 
gram of the overall system is shown in figure 4. 

Analog, synchro, frequency, and ARINC 561 format 
data signals received from the TAMMS components are 
sent, via signal conditioning and interface electronics, 
to the PCM system which uses a remote multiplexer/ 
demultiplexer unit (RMDU) for processing data. Synchro 
and analog data are recorded at 1 6- and 12-bit accuracies, 
respectively. Selected analog signals are filtered by using 
8 pole, 6 zero elliptic low-pass filters. Signal cutoff fre- 
quencies varied accordingly with the response character- 
istics of the respective instrumentation (e.g., EG&G 137 
hygrometer, 2 Hz; total air temperature sensor, 5 Hz; ver- 
tical velocity and fast-response measurements of chemi- 
cal species, 10 Hz). Data are recorded at 66 Hz by a 
dedicated 386 personal computer onto a dual disk drive. 
The dual drive allows a full data disk to be swapped 
without loss of data. A flight data recorder is also used as 
a storage backup. A second 386 personal computer per- 
forms system control, data analysis, and display tasks in 
real time. 

4. Analysis 

4.1. Static Pressure Error 

The flow field in front of an airplane is influenced 
for a relatively large distance at subsonic speeds. Far 
upstream from the airplane, where fuselage effects on the 
streamlines are minimal, the total pressure is the sum of 
the free-stream static and dynamic pressures, 

Pt = P + V c (?) 

The shape of the fuselage nose, primarily, determines the 
distance ahead of the airplane that the flow field is 
altered. Because of the blunt nose of the Electra airplane, 
the streamlines of the oncoming air must diverge rapidly; 
thus, differences are created between the measured and 
free-stream pressures. Total pressure p t is conserved for 
any point along the aircraft body for subsonic flight as 
given by 

P t = P + <4 ( 8 ) 

where the prime denotes local values. The static pressure 
error A p defined as 

Ap = p ~ p (9) 

is substituted into equation (8) to get 

p t = p + Ap + q c (10) 
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Subtracting equation (7) from equation (10) gives 

q c = <i c ~^P 


(ii) 


during flight, angle of attack is usually correlated with 
nondimensional position error: 


The magnitude of Ap can be substantial for very 
blunt fuselages. Wind tunnel tests have previously dem- 
onstrated that the magnitude of static pressure error, as a 
ratio of Apfq c , is influenced by nose shape and distance 
upstream of the nose (Letko 1 947). The presence of the 
airplane fuselage, therefore, causes an increase in the 
local static pressure as well as a decrease in the local 
dynamic pressure. Changes in the local pressure field 
may also result from airflow distortion induced by the 
presence of the pitot-static tube. The determination of the 
overall static pressure error is a requisite and can be 
determined through flight calibration tests (Gracey 1980; 
Lenschow 1986; Dommasch, Sherby, and Connolly 
1961). 

An additional source of error in pitot-static systems 
arises from resonance in the lines between the pressure 
port and the pressure transducer (Iberall 1950) and, 
therefore, is a function of the length and diameter of tub- 
ing and air density (Nicholls 1982). In the NASA Electra 
configuration, the length and diameter of tubing are 
approximately 0.75 m and 0.25 cm, respectively; thus, 
the resonant frequency is ensured to be above 1 0 Hz. 

Gracey (1958 and 1980), Etkin (1959), and 
Dommasch, Sherby, and Connolly (1961) primarily 
associate changes in airflow, which influences Aplq CJ 
with Mach number, angle of attack, and/or lift coeffi- 
cient. The Mach number is given as 


M 


2 


2 

Y-l 



(12) 


where q c and p are, respectively, the dynamic pressure 
and free-stream static pressure (defined earlier). The 
ratio of specific heats y, which has a slight dependence 
on water vapor (List 1971), is given by 


y = y o (l- 0.07 \29-^ (13) 

where e is partial vapor pressure and y 0 = 1 .4028 for dry 
air. 


As mentioned, Aplq c is influenced by the lift coeffi- 
cient C^, which is defined as 


C 


L ~ 


L 

?c S 


(14) 


^ = /(a) (15) 

<lc 

This correlation is possible because C L and angle of 
attack are linearly related for small angles and a given 
airplane configuration (Gracey 1980); this allows angle 
of attack to be substituted for C L in these conditions. 
Thus, from equation (14) the nondimensional static pres- 
sure error at the static port varies with angle of attack. 
Altering the airplane configuration also affects the flow 
field and position error. Pressure errors due to accelera- 
tion are relatively small and become negligible with 
flight leg scenarios which are primarily kept at steady 
conditions and below M = 0.5. Consideration of a much 
wider dynamic range of flight parameters is recognized 
for the calibration of high performance aircraft and 
expanded range of flight conditions (Haering 1990; 
Bjarke and Ehernberger 1989). 

4.2. Static Air Temperature 

The static air temperature (i.e., the temperature of 
the static, undisturbed air through which the airplane is 
traveling) is actually the last of four temperature values 
that must be determined for any given point. The four 
temperatures in order of derivation are (1) measured 
temperature T m , (2) recovery temperature T n (3) total 
temperature T t , and (4) static temperature T s . These 
parameters are discussed in detail by Stickney, Shedlov, 
and Thompson (1990). For our purposes a short sum- 
mary will suffice. 

Static air temperature T s is related to the total tem- 
perature T t by 


t > = t { x+T T m 2 i) (16) 

where M l denotes the local Mach number as determined 
by local values of static and dynamic pressure. Estima- 
tions of the difference in local static and dynamic pres- 
sures were made for the Rosemount total temperature 
probe located 0.6 m downstream from the pressure ports. 
The increase of the local static pressure at the Rosemount 
sensor was thus estimated to be approximately 1 0 percent 
of the free-stream dynamic pressure. 

The recovery temperature T r is related to T t by 


where L is lift which is equivalent to the weight of air- 
plane for level flight and S is wing area. Because of the 
difficulty in determining and monitoring aircraft weight 


T < = 


1 -T) 


(17) 
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where T| is the recovery correction given by 


(l-r)[(Y-l)/2]Mf 

T1 = 08) 

1 +[(y-\)/2]M, 

where r is the recovery factor for the sensor. The recov- 
ery factor is a function of the aerodynamics of the sensor 
and its location on the aircraft and must therefore be 
determined experimentally. For each data point as the 
aircraft passed the tower, r was determined from 



where T s ref is the reference static temperature as deter- 
mined from tower measurements interpolated to aircraft 
altitude. Finally the measured temperature T m is related 
to T r by 


T = T -T 

r m x 


tion error and T s for the TAMMS platform. A description 
of the test setup and procedure is in section 5.1. 

In-flight maneuvers performed remotely from the 
tower are also included in the calibration procedure. 
These maneuvers modulate var ious airflow configura- 
tions relative to the airplane and are required to correct 
flight measurements and evaluate the overall perfor- 
mance of TAMMS. Procedures call for these maneuvers 
to be performed for the initial calibration and determina- 
tion of systematic system errors. The need to account for 
any changes in flow vane characteristics and ensure data 
integrity also requires that these maneuvers be performed 
on a regular basis on data flights. The various in-flight 
maneuvers utilized for air motion system calibrations are 
discussed by Axford (1968), Telford and Wagner (1974), 
Lenschow et al. (1978), and Lenschow (1986). The pro- 
cedures applied to the overall calibration of TAMMS 
consist of speed variation, skid, pitch, and reversed head- 
ing maneuvers, and are discussed in sections 5.2, 5.3, 
6.1, and 6.2. 


where T x is an error due to self-heating effects and is 
generally negligible. True airspeed, as indicated in equa- 
tion (3), contributes in the calculation of all wind vector 
components. 

True airspeed U a is expressed as 

u 2 a = y RT S M 2 (21) 

where R is gas constant. It is apparent from equations (3), 
(9), (11), (12), and (21) that the accuracy of air motion 
measurements is dependent on the accurate determina- 
tion of T s and free- stream impact and static pressures. 
The temperature and pressure measurements used in the 
calculation of U a , however, must be determined from 
measurements that are characteristic of the free-stream or 
undisturbed air. Since these measurements were obtained 
from sensors located on the gust probe, they must be cor- 
rected for kinetic energy effects and position error due to 
the disturbance of the airplane body to the free-stream 
airflow. These corrections for the TAMMS are deter- 
mined through dedicated in-flight maneuvers and valida- 
tion tests which are discussed in detail in sections 5 
and 6. 

5. Flight Calibration Techniques 

Techniques employed to calibrate sensors and deter- 
mine systematic errors encountered by air motion mea- 
surement systems are well known and documented 
(Shrager 1964; DeLeo and Hagen 1966; Anon. 1971; 
Gracey 1980; DeAnda 1981; and Brown 1985). A varia- 
tion of the tower flyby method with radar tracking was 
implemented for the calibration of static pressure posi- 


5.1. Tower Flyby 

A 91.5-m tower located at the Wallops Flight Facil- 
ity was used for flight calibration of the TAMMS plat- 
form (fig. 5). A tracking C-band radar located 3483 m 
from the instrumented tower was used to record precise 
three-dimensional position data of the Electra airplane 
during the TAMMS calibration flight. The radar employs 
an rf amplifier and is designed specifically for precision, 
long-range tracking. The Electra was outfitted with a 
transponder for optimal precision tracking. The range 
and angle precision of the radar is 2.75 m (rms) and 
0.05 mils (rms), respectively. This translates to approxi- 
mately ±0.2 m in aircraft altitude. The radar transmit- 
ting and receiving frequency range is 5400 to 5900 MHz 
with a minimum range of 550 m. 

The tower was instrumented to measure absolute 
pressure, air temperature, and dew point. Air temperature 
measurements were made by thermistor sensors with an 
accuracy of ±0.2°C at 91.5 m, 76.2 m, and ground lev- 
els. Static air pressure and dew point measurements were 
made at 91.5 m. The absolute atmospheric pressure sen- 
sor consisted of a Paroscientific pressure transducer (0 to 
1035 hPa) housed in a thermally controlled enclosure for 
environmental stability and is identical to the transducer 
used on the TAMMS platform. (See section 3.) Preced- 
ing and following the calibration flight, further correla- 
tion between the airplane and tower pressure sensors was 
obtained by direct comparison to a SONIX digital pres- 
sure system developed by Pressure Systems, Inc., to meet 
calibration standards. 

The airplane path was straight, level, unaccelerated, 
and just below the top of the tower for approximately 
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30 s before and after passage of the tower. Tower flyby 
calibration data were taken during 21 passes at 5 differ- 
ent speeds. A clean, aerodynamic profile was maintained 
(i.e., no flaps) at all times. Passes were made with the 
flight path offset approximately 2 to 3 wingspans 
(100 m) from the tower. Data taken on the airplane were 
averaged at 10 Hz to match the radar sampling rate. Ref- 
erence values for air temperature T s re f and pressure p S T& f 
were determined for comparison with values obtained 
from the TAMMS platform within 3 s of passage of the 
tower. These reference values were based on the height 
data obtained from the radar and the temperature and 
pressure reference measurements made on the tower. 
Because pitch attitude changed with airspeed of the air- 
plane during the tower flybys, the reference height 
(determined from radar tracking) was geometrically cor- 
rected for the difference between the vertical position of 
the radar transponder and static pressure sensor. 

Letko (1947) published wind tunnel data indicating 
the variation of A p/q c with x ID for a body with a blunt- 
shaped nose (jc is distance of static port in front of nose 
and D is fuselage diameter). Wind tunnel data shown in 
figure 6 indicate A plq c is constant during level flight for 
Mach numbers within the Electra flight envelope and xID 
ratio (O’Bryan, Danforth, and Johnston 1955). Roe 
(1951), however, also presents data (fig. 7) that illustrate 
a slight Ap/q c relationship to angle of attack; thus, tower 
flyby data for the TAMMS configuration (fig. 8) support 
a relationship between Aplq c and angle of attack. 

Static pressure error (difference between TAMMS 
static pressure and free- stream static pressure interpo- 
lated from tower to airplane height), obtained from tower 
flybys, ranged from approximately 3 hPa at Mach num- 
ber 0.2 to 9 hPa at Mach number 0.45. An approximate 
value of 0.075 was determined for A p!q c with a slight 
correlation to angle of attack as indicated by the tower 
flyby data shown in figure 8. These results are highly 
correlative to the wind tunnel data in figures 6 and 7. 

Figure 9 shows the experimentally obtained values 
of the recovery factor for the temperature probe as a 
function of M /. These data indicate that although a 
second-order curve is shown to describe the variation of r 
with airspeed, a value of 0.98 could be used over the 
speed range of the airplane with minimal error. Calcu- 
lated total and static air temperatures as a function of A// 
are shown in figure 10. The obvious correlation of air- 
speed with T t is not present, as should be expected, for 
the calculated T s data curve. The difference between the 
reference tower temperature and the airborne static tem- 
perature measurements for each pass of the tower was 
determined and indicated no dependence on airspeed. 
Although no speed dependence is noticeable, the air- 
borne measurements were found to be ~0.3°C lower than 


the reference values, which are near measurement accu- 
racies of the instrumentation. 

5.2. Speed Variation Maneuver 

The speed variation maneuver was the first of sev- 
eral in-flight maneuvers required to correct measure- 
ments of flow angles due to the distortion airflow around 
the airplane. The maneuver consists of flying at a con- 
stant heading and altitude while varying airspeed through 
the speed range of the airplane. This maneuver modu- 
lates angle of attack a and pitch angle 0 due to the rela- 
tion of lift to a and airspeed and allows the alpha vane to 
be directly calibrated to 0, which is obtained from the 
INS. Because of upwash effects in front of the airplane, a 
correction or sensitivity factor is found to be 0.75 for a 
as indicated by the slope of the plot of a versus 0 shown 
in figure 1 1 . The sensitivity of the alpha vane is not only 
affected by upwash, but a zero offset is also determined 
to be 1.35° from the same data set. This offset includes 
any rotational misalignment of the alpha vane synchro 
repeater caused by the nose boom configuration and 
installation. This procedure was repeated during mis- 
sions to compensate for shifts in offset and sensitivity of 
the flow angle vanes. 

These maneuvers are also used to check the depen- 
dency of T s on airspeed. Adiabatic effects on the temper- 
ature sensor (i.e., measured temperature) vary with 
airspeed. The removal of these effects is done through 
computations which require a known recovery factor for 
the sensor. Errors in determining this value or the mea- 
surement of airspeed show a periodic relationship 
between T s and airspeed; this indicates a dependence of 
T s on airspeed. The time series plot of T s and airspeed 
shown in figure 12, however, indicates no such relation- 
ship and reinforces the validity of airspeed and tempera- 
ture sensor calibrations. 

5.3. Skid Maneuver (Yaw or Sideslip) 

Airplane heading \| /, horizontal velocity components, 
and sideslip (J are modulated during the skid maneuver 
while altitude and roll are kept constant. Skid maneuvers 
are necessary for in-flight calibration of the sensitivity 
factor and estimation of errors associated with the beta 
vane. Sensitivity and offset values are unique to each 
flow vane installation and must be determined accord- 
ingly. Sideslip flight calibrations, however, cannot 
assume that the horizontal winds are negligible in con- 
trast to the calibration of the angle of attack when vertical 
winds are assumed to be zero. Tjemstrom and Friehe 
(1991) describe a flight calibration technique employed 
for sideslip, which assumes for homogeneous low- 
turbulence conditions that all variations in the computed 
horizontal wind are induced by the skid maneuver only. 
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The sensitivity factor determined under homogeneous 
and low-turbulence conditions during skid maneuvers is 
assumed to be highly representative of the beta vane. The 
offset value, which may include a bias error related to \|/ 
from the INS alignment, however, is determined by 
adjusting the offset value until the cross-track wind com- 
ponent is equivalent in magnitude for both directions of a 
reversed heading maneuver. Calibration errors associated 
with the beta vane are indicated by a periodic variation in 
computed horizontal velocity during skid maneuvers. 
The beta vane sensitivity factor was found to be 0.81 
with an offset of 1.10°. These values rendered horizontal 
wind field measurements of less than 10 percent of 
U a ( sin p), which meet the criteria for satisfactory system 
performance. 

Sideslip, similar to angle of attack, influences static 
pressure measurements. A time series plot of the aircraft 
sideslip motion, measured static pressure, and corrected 
static pressure for a skid maneuver is shown in figure 13. 
The large sinusoidal curve indicates the sideslip angle as 
measured by the rotating vane. A split scale is used, for 
data presentation purposes, to offset the measured static 
pressure ( p\ upper curve) from the final static pressure 
values (p, lower curve) which are determined through 
calibration as a function of P and q c . 

Data from sideslip maneuvers indicate an asymmet- 
ric relationship of sideslip with calculated static air tem- 
perature T' s . A time history plot of sideslip, T', and 
corrected T s is presented (in similar fashion as the previ- 
ous static pressure data set) in figure 14. The variation in 
T' s was removed for positive sideslip angles by using the 
corrected static pressure to determine M l values. This 
approach, however, was not sufficient to remove the 
larger variations in T' during negative sideslip angles. 
The nonsymmetrical relation of T' with sideslip is 
attributed to the nose boom configuration (fig. 3), which 
shows the total temperature probe mounted horizontally 
on the starboard side of the boom. This position produces 
a shielding effect during a negative sideslip angle. An 
empirical multiplication factor was determined and 
applied to previously determined static pressure correc- 
tions to yield local static pressure conditions and subse- 
quently Mi encountered by the total temperature probe 
during negative sideslip. The resultant T s with the major- 
ity of the asymmetric effects from shielding removed is 
indicated by the lower curve in figure 14. 

6. Flight Validation Maneuvers 

6.1. Pitch Maneuver 

Vertical airplane velocity, airspeed, and angle of 
attack are modulated during a pitch maneuver, whereas 
roll and sideslip are kept constant. The purpose of the 


pitch maneuver is to indicate any errors in the determina- 
tion of vertical winds w or, conversely, demonstrate the 
overall ability of the air motion system to remove vertical 
airplane motion and make ambient vertical wind mea- 
surements. Lenschow (1986) suggests a period of 10 to 
20 s for the pitching motion with an airplane vertical 
velocity of 2.5 to 4 m-s” 1 for this maneuver. The airplane 
vertical velocity and ambient vertical wind determined 
from a pitch maneuver are plotted for comparison in fig- 
ure 15. The amplitude of airplane vertical velocity oscil- 
lation is varied between 6 and 20 m-s' 1 with a 20- to 30-s 
period. The general performance criteria for an air 
motion system are that errors in calculated vertical winds 
be less than 10 percent of the airplane vertical velocity. 
These criteria, as illustrated, are easily met during the test 
conditions. 

6.2. Reversed Heading Maneuver 

The reversed heading maneuver demonstrates the 
ability of the air motion system to measure the ambient 
horizontal wind components u and v by modulating air- 
speed and sideslip angle errors. From equations (3), \j /, 
u p , and v p are obtained through the INS leaving U a and [} 
as determined by the pitot-static tube, T t sensor, and 
beta vane measurement; thus, the calibrations of these 
sensors are critical and directly influence horizontal wind 
measurements. 

Reversed headings are performed by flying at a con- 
stant altitude and heading for approximately 3 min and 
then reversing the course in a manner in which the return 
course passes through the same air mass. These maneu- 
vers are performed above the boundary layer in a homo- 
geneous air mass to ensure the sampling of the same 
meteorological conditions as realistically possible. A 
time history plot of along-track winds during a reversed 
heading maneuver performed along the prevailing wind 
direction is shown in figure 16(a) and demonstrates the 
ability of the system to make accurate, true airspeed mea- 
surements. A reversed heading is immediately performed 
adjacent to the first and perpendicular to the prevailing 
wind. Cross-track winds obtained in these conditions are 
shown in figure 16(b) and are similar, as expected, to the 
along-track winds. This similarity indicates that the sum 
of sideslip and heading has been adequately determined. 
Horizontal wind, obtained by vector addition of along- 
and cross-track winds from a previous flight track is 
shown in figure 16(c). 

7. Concluding Remarks 

The calibration of a turbulent air motion system from 
in-flight and tower flyby maneuvers using a tracking 
C-band radar has been discussed. A description of the 
gust probe measurement technique and instrument 
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package, flight calibration procedures, and results were 
included. The primary goals that have been addressed 
consist of: (1) defining any dependence of static air tem- 
perature to airspeed, (2) determining systematic static 
pressure position error as related to platform configura- 
tion, and (3) demonstrating the overall system consis- 
tency to make accurate three-dimensional ambient wind 
field measurements. 

The total air temperature probe recovery factor of 
0.98 was determined experimentally through analysis of 
the tower flyby calibration data. Comparison of airplane- 
derived static air temperature T s with tower air tempera- 
ture data and speed variation maneuvers performed 
remotely from the tower illustrated the lack of depen- 
dence between T s and airspeed; consequently, the 
calibration of the air temperature probe was verified. 
Further flight tests (skid maneuvers), however, indicated 
a shadowing effect on the total air temperature probe 
induced by the pitot-static tube mounting configuration 
on the nose boom. This shadow effect was successfully 
removed through the application of a numerical filter to 
the air temperature data. 

Calibration data acquired for the determination of 
static pressure position error of the gust probe were 
found to closely approximate data obtained from previ- 
ously documented wind tunnel tests. Analysis of tower 


and flight data showed a high correlation of position 
error to impact pressure with a slight dependence on 
angle of attack and/or lift coefficient C L . Because of the 
difficulty in obtaining a reasonable approximation of air- 
plane weight during flight and the linear relationship of 
angle of attack with C L for small angles, angle of attack 
was substituted for C L . As a result, static pressure error 
was determined to be defined adequately as a function of 
impact pressure and angle of attack. Pitch, skid, and 
reversed heading maneuvers demonstrated the overall 
ability of the air motion system to remove aircraft motion 
and make ambient vertical and horizontal wind measure- 
ments. The general rule applied to airborne platforms 
used to make air motion measurements states that the cal- 
culated vertical winds should be <10 percent of the plat- 
form vertical velocity and calculated horizontal winds 
should be <10 percent of U a ( sin P) (where U a is true air- 
speed and P is sideslip angle). These criteria have been 
met and provide ample confidence in the system to 
provide reliable measurements of the ambient three- 
component wind field during straight and level flight 
conditions. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
October 4, 1996 
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Appendix 


Specifications of Lockheed 188 Electra Airplane 

The Lockheed 188 Electra airplane based at Wallops Flight Facility and designated NASA 429 is a low-wing, 
medium altitude, moderate range transport airplane powered by four Allison turboprop engines.The specifications of the 
airplane are given as follows: 


Wing: 

Area, m 2 120.8 

Span, m 30.2 

Fuselage: 

Diameter, m 3 45 

Overall length (excluding boom and pitot-static tube), m 32.25 

Weight: 

Takeoff gross weight, kg 52 730 

Payload capacity, kg 10450 

Maximum fuel loading, kg 15 450 

Performance: 

True airspeed, m-s' 1 90-160 

Flight endurance at 155 m-s -1 , hr 5.5 

Maximum operating altitude, km 7 

Cruising altitude, km 4.5-7 


Altitude, km 

Endurance, hr 

Range, km 

Airspeed, m-s 1 

0.15-3 

6 

2600 

125 

3-6 

7 

2800 

125 

6-7 

7 

3100 

130 
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Table I. TAMMS Sensor Characteristics 


Parameters 

Sensor 

Accuracy 

Resolution 

Response 

Dew point 

EG&G 137 hygrometer 

±0.5°C 

0.1 °c 

2-10 sec 

Dew point 

AIR-LA- 1 AC Lyman- Alpha 
hygrometer 

±0.6°C 

0.03 °C 

2 ms 

Air temperature 

Rosemount model 102 non-deiced 
series sensor 

±0.2°C 

0.006 °C 

a 2 Hz 

Flow angles 

Free-rotating balsa vanes (synchro 
resolver) 

±0.2 °C 

0.005 °C 

25-50 Hz 

Static pressure 

Paroscientific Digiquartz absolute 
pressure transducer 

±0.21 hPa 

0.052 hPa 

1 ms 

Differential pressure 

Paroscientific Digiquartz differential 
pressure transducer 

±0.083 hPa 

0.021 hPa 

1 ms 

Pressure altitude 

Rosemount 124 IB capacitive 
transducer 

±0.4 percent 
(0 to 4.5 km) 

1 m at sea level 

15 ms 


a Extended to =8 Hz with numerical filter (Ritter et al. 1987). 
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Front view 


Figure 1 . Coordinate system, airplane attitude angles, and reference axes used to calculate air velocity components. 
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WI-84-368-6 

Figure 2. Lockheed 1 88 Electra airplane (NASA 429) at Wallops Flight Facility. 



Figure 3. TAMMS gust probe. 
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Figure 4. Block diagram of TAMMS data acquisition system. 
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Mach number 


Figure 6. Wind tunnel data illustrating variation of static pressure error as influenced by Mach number at different 
distances ahead of aircraft fuselage (O’Bryan, Danforth, and Johnston 1955). 



o .4 .8 1.2 1.6 2.0 


Figure 7. Wind tunnel data showing effect of angle of attack on static pressure error forward of airplane with blunt- 
shaped nose similar to Electra airplane. (Data from Roe 1951.) 




Figure 8. Static pressure error determined from tower flyby test data as function of angle of attack. 
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Figure 10. Tower flyby data indicating dependent nature of total air temperature and required isolation of static air 
temperature on Mach number. 
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Static pressure, 
p\ mb 


Static pressure, 
/?, mb 


Figure 13. Time series plot of sideslip, measured static pressure p\ and corrected static pressure p obtained during skid 
maneuver. 



Figure 14. Time series plot of sideslip, calculated static air temperature T' s , and corrected static air temperature T 
obtained during skid maneuver. 
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Figure 15. Time series plot of aircraft vertical velocity and vertical wind component as determined by TAMMS during 
pitch maneuver. 
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(a) Along-track wind. 


(b) Cross-track wind. 
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(c) Horizontal wind. 

Figure 16. Plots of winds as determined by vector addition. 
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